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Abstract 

A family of non-equilibrium statistical operators (NSO) is introduced 
which differ by the system lifetime distribution over which the quasi- 
equilibrium (relevant) distribution is averaged. This changes the form of 
the source in the Liouville equation, as well as the expressions for the 
kinetic coefficients, average fluxes, and kinetic equations obtained with 
use of NSO. It is possible to choose a class of lifetime distributions for 
which thermodynamic limiting transition and to tend to infinity of av- 
erage lifetime of system is reduced to the result received at exponential 
distribution for lifetime, used by Zubarev. However there is also other ex- 
tensive class of realistic distributions of lifetime of system for which and 
after to approach to infinity of average lifetime of system non-equilibrium 
properties essentially change. For some distributions the effect of "finite 
memory" when only the limited interval of the past influence on be- 
haviour of system is observed. It is shown, how it is possible to spend 
specification the description of effects of memory within the limits of 
NSO method, more detailed account of influence on evolution of system 
of quickly varying variables through the specified and expanded form of 
density of function of distribution of lifetime. The account of character 
of history of the system, features of its conduct in the past, can have sub- 
stantial influence on non-equilibrium conduct of the system in a present 
moment time. 
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1 Introduction 

One of the most fruitful and successful ways of development of the description of 
the non-equilibrium phenomena are served by a method of the non-equilibrium 
statistical operator (NSO) [BO- In work [3] new interpretation of a method 
of the NSO is given, in which operation of taking of invariant part [TJ [5] or 
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use auxiliary "weight function" (in terminology [H [5]) in NSO are treated 
as averaging of quasi-equilibrium statistical operator on distribution of past 
lifetime of system. This approach adjust with the operations spent in the 
general theory of random processes, in the renewal theory, and also with the 
lead Zubarev in work as [2] reception NSO by means of averaging on the initial 
moment of time. 

This treatment of NSO gives to the procedure looking before formal, physi- 
cal sense of the account of causality and allocation of a real finite time interval 
in which there is a given physical system. New interpretation leads to vari- 
ous directions of development of NSO method which is compared, for example, 
with Prigogine's [5] approach, introduction of the operator of internal time, 
irreversibility at microscopical level. 

In Kirkwood's works [7] it was noticed, that the system state in time present 
situation depends on all previous evolution of the non-equilibrium processes 
developing it. For example, in real crystals it is held in remembrance their 
formation in various sorts "defects" (dispositions etc.), reflected in structure 
of the crystals. Changing conditions of formation of crystals, we can change 
their properties and create new materials. In works [H [5] it is specified, that 
it is possible to use many " weight functions" . Any form of density of lifetime 
distribution gives a chance to write down a source of general view in dynamic 
Liouville equation which thus becomes, specified Boltzmann and Prigogine [H 
03 [5] , and contains dissipative items. 

If in Zubarev's works CD H] the linear form of a source corresponding limit- 
ing exponential distribution for lifetime is used other expressions for density of 
lifetime distribution give fuller and exact analogues of " integrals of collisions" . 
The obvious account of violation of time symmetry (through finiteness of life- 
time, the beginning, the end and irreversibility of a life) is entered. Besides 
communication with the theory of queues, reliability theory, the management 
theory, the information theory etc., in offered work the physical consequences 
connected with fundamental physical problems are reflected. 

The formalism follows from the physical matter, for example, from finite- 
ness of lifetime of real physical systems (it is possible to result many examples 
of problems in which it is necessary to consider systems of the finite sizes with 
finite lifetime). Generally the description of non-equilibrium systems repre- 
sents the self-coordinated problem: definition of lifetime through interaction 
of system with environment CQ , dynamics of the operators characterizing non- 
equilibrium processes, and substitution of found average lifetime in NSO, defini- 
tion of non-equilibrium physical characteristics, depending on system lifetime. 

In work [8] irreversible transfer equations are received in assumptions of 
coarsening of the distributions, a certain choice of macroscopical variables and 
the analysis of division of time scales of the description (last circumstance was 
marked in [3]). Importance and necessity of the analysis of the time scales 
playing a fundamental role in the description of macroscopical dynamics of 
system is underlined. Evolution of slow degrees of freedom is described by 
Markovian equations. Thus the time scale on which observable variables evolve, 
should be much more time of memory on which the residual effects brought by 
irrelevant degrees of freedom are considered. 
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Otherwise effects of memory play an essential role. Memory time is esti- 
mated in work [8 for Boltzmann equation. In the present work the consider- 
ation subject is made by situations when it is necessary to consider effects of 
memory. Examples of such situations are given in [S]. We will notice, that 
the projective methods used in [5], do not consider distribution on lifetime of 
system (that is noted in [5]) on which method NSO is based. 

In work [10] it is shown, in what consequences for non-equilibrium properties 
of system results change of lifetime distribution of system for systems of the 
limited volume with finite lifetime. In the present work are considered also 
infinitely greater systems with infinite average lifetime. 

2 New interpretation of NSO 

In [3] the Nonequilibrium Statistical Operator introduced by Zubarev [U [2] 
rewritten as 



p oo 

lng{t) = / p q (u)lng q (t — u, —u)du, lng q (t,Q) = — $(t) — ~S^F n {t)P n \ 
Jo 

Ing^tM) = e^ H / a hng q (t,0)e^ H /^; *(t) - lnSpexp^i^i^}, 

71 

where H is hamiltonian, lng(t) is the logarithm of the NSO in Zubarev's form, 
lng q (t,0) is the logarithm of the quasi-equilibrium (or relevant); the first time 
argument indicates the time dependence of the values of the thermodynamic 
parameters F m ; the second time argument ti in Q q {t\,t<i) denotes the time 
dependence through the Heizenberg representation for dynamical variables P m 
from which g q (t, 0) can depend [TJ [U [31 [H [5] . In [3] the auxiliary weight 
function p q (u) = eexp{— eu} was interpreted as the probability distribution 
of lifetime density of a system. T is random variables of lifetime from the 
moment to of its birth till the current moment t; e _1 = (t — to); (t — to) = (T), 
where (r) = J up q (u)du is average lifetime of the system. This time period 
can be called the time period of getting information about system from its 
past. Instead of the exponential distribution p q {u) in (1) any other sample 
distribution could be taken. This fact was marked in [3] and [H [5] (where the 
distribution density p q (u) is called auxiliary weight function w(t, t')). From the 
complete group of solutions of Liouville equation (symmetric in time) the subset 
of retarded "unilateral" in time solutions is selected by means of introducing a 
source in the Liouville equation 

+ iLg(t) = -e(g(t) - g q (t, 0)) = J, 

which tends to zero (value e — + 0) after thermodynamic limiting transition. 
Here L is Liouville operator; iL = -{H,g} = £ fe [f^ - §if& H is 
Hamilton function, pk and qk are pulses and coordinates of particles; {...} is 
Poisson bracket. In [11] it was noted that the role of the form of the source term 
in the Liouville equation in NSO method has never been investigated. In [12] it 
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is stated that the exponential distribution is the only one which possesses the 
Markovian property of the absence of contagion, that is whatever is the actual 
age of a system, the remaining time does not depend on the past and has the 
same distribution as the lifetime itself. It is known [TJ [2J [31 HI [5] that the Liou- 
ville equation for NSO contains the source J — J zu b — —e[lng(t) — lng q (t, 0)] 
which becomes vanishingly small after taking the thermodynamic limit and 
setting e —> 0, which in the spirit of the paper [T| corresponds to the infinitely 
large lifetime value of an infinitely large system. For a system with finite size 
this source is not equal to zero. In [5] this term enters the modified Liou- 
ville operator and coincides with the form of Liouville equation suggested by 
Prigogine [B] (the Boltzmann-Prigogine symmetry), when the irreversibility is 
entered in the theory on the microscopic level. We note that the form of NSO 
by Zubarev cast in [3] corresponds to the main idea of [6] in which one sets 
to the distribution function g (g q in Zubarev's approach) which evolves ac- 
cording to the classical mechanics laws, the coarse distribution function g (g(t) 
in the case of Zubarev's NSO) whose evolution is described probabilistically 
since one perform an averaging with the probability density p q (u). The same 
approach (but instead of the time averaging the spatial averaging was taken) 
was performed in [13) . 

Besides the Zubarev's form of NSO [HH], NSO Green-Mori form [H US] 
is known, where one assumes the auxiliary weight function [3] to be equal 
W(t,f) = 1 - (t - t')/T;w(t,t') = dW(t,f)/dt' = l/r;r = t - t . After 
averaging one sets r — > oo. This situation at p q (u = t — to) = w(t,t' = to) 
coincides with the uniform lifetime distribution. The source in the Liouville 
equation takes the form J — lng q /r. In [1] this form of NSO is compared to 
the Zubarev's form. 

One could name many (no less than 1000) examples of explicit defining of 
the function p q (u). Every definition implies some specific form of the source 
term J in the Liouville equation, some specific form of the modified Liouville 
operator and NSO. Thus the family of NSO is defined. If distribution p q (u) 
contains n parameters, it is possible to write down n equations for their expres- 
sion through the parameters of the system. From other side, they are expressed 
through the moments of lifetime. There is the problem of optimum choice of 
function p q (u) and NSO. 

3 Modifications to the nonequilibrium descrip- 
tion 

Let's consider now, what consequences follow from such interpretation of NSO. 
3.1 Families of NSO 

Setting various distributions for past lifetime of the system, we receive a way of 
recording of families of NSO. Class of NSO from this family will be connected 
with a class of distributions for lifetime (taken, for example, from the stochastic 
theory of storage processes, the theory of queues etc.) and with relaxation 
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properties of that class of physical systems which is investigated. The general 
expression for NSO with any distribution 

/•OO 

lng(t) = / p q (u)lnQq(t — u, —u)du = (1) 
Jo 

i n \ f°° i f r \ ^ din g(t — u, —u) 
= lnQ q {t,0)-J (J p q (u)du) — du, 

where integration by parts in time is carried out at J p q {y)dy\ y=0 = — 1; 
J ' Pq(y)dy\ y ^ oc> = 0; at p q (y) = eexp{-ey};e = l/(T), the expression (1) 
passes in NSO from [TJ E]- In [12] it is shown, how from random process 
X(t), corresponding to evolution of quasi-equilibrium system, it is possible to 
construct set of new processes, introducing the randomized operational time. 
It is supposed, that to each value t > there corresponds a random value 
T(t) with the distribution p q (y). The new stochastic kernel of distribution 
of a random variable X(T(t)) is defined by equality of a kind (1). Random 
variables X(T(t)) form new random process which, generally speaking, need not 
to be of Markovian type any more. Each moment of time t of "frozen" quasi- 
equilibrium system is considered as a random variable T(t) the termination of 
lifetime with distribution p q {y). Any moment of lifetime can be with certain 
probability the last. That the interval t — ta = y was enough large (that became 
insignificant details of an initial condition as dependence on the initial moment 
to is nonphysical p][2]), it is possible to introduce the minimal lifetime T m i n = 
T% and to integrate in (fl]) on an interval (Ti,oo). It results to the change 
of the normalization density of distribution p q (y). For example, the function 
P q {y) = eexp{-ey} will be replaced by p q (y) = e exp{eri-ey}, y > Ti;p q (y) = 
Q,y < Ti. The under limit of integration in {l} by Ti — > is equal 0. It is 
possible to choose p q (y) = Cf(y),y < h\p q {y) = eexp{~ey},y >h;C = (1 - 
exp{— sti})/(J Q 1 f(y)dy). The function f{y) can be taken from models of the 
theory of queues, the stochastic theory of storage and other sources estimating 
the lifetime distribution for small times (for example [H [H EEE1 123 ) ■ The 
value t\ can be found from results of work [18] . It is possible to specify many 
concrete expressions for lifetime distribution of system, each of which possesses 
own advantages. To each of these expressions there corresponds own form of 
a source in Liouville equation for the nonequilibrium statistical operator. In 
general case any functions p q (u) the source is: 

J = Pq (0) In Qq (t, 0) + £ ^M(lne q (t ~ y, -y))dy (2) 

(when values J? g (0) disperse, it is necessary to choose the under limit of inte- 
gration equal not to zero, and T m i n ). Such approach corresponds to the form 
of dynamic Liouville equation in the form of Boltzmann-Bogoliubov-Prigogine 
[HO [6], containing dissipative items. 

Thus the operations of taking of invariant part pQ, averaging on initial 
conditions [2], temporary coarse-graining [7], choose of the direction of time 
[H[S], are replaced by averaging on lifetime distribution. 
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The physical sense of averaging on introduced lifetime distribution of quasi- 
cquilibrium system as it was already marked, consists in the obvious account 
of infringement of time symmetry and loss (reduction accessible) the informa- 
tion connected with this infringement, that is shown in occurrence the value 
of average of entropy production (AS(t)) not equal to zero, obviously reflect- 
ing fluctuation-dissipative processes at the real irreversible phenomena in non- 
equilibrium systems. The correlations received in the present section general- 
ize results of statistical non-equilibrium thermodynamics [T] [5] and informa- 
tion statistical thermodynamics [H [S] as instead of weight function of a form 
e exp{st'} contain density of probability of lifetime of quasi-cquilibrium system 
which as it was already marked, can not coincide with exponential distribution 
(in the latter case it coincides with weight function from [HI!])- For example, 
for system with n classes of ergodic states limiting exponential distribution is 
replaced with the general Erlang. In research of lifetimes of complex systems 
it is possible to involve many results of the theory of reliability, the theory of 
queues, the stochastic theory of storage processes, theory of Markov renewal, 
the theory of semi-Markov processes etc. 

It is essentially that e ^ 0. The thermodynamic limiting transition is not 
performed, and actually important for many physical phenomena dependence 
on the size of system are considered. We assume e and (r) to be finite values. 
Thus the Liouville equation for g(t) contains a finite source. The assumption 
about fmitcness of lifetime breaks temporary symmetry. And such approach 
(introduction p q (y), averaging on it) can be considered as completing the de- 
scription of works [U [2] • 

In work [18] lifetimes of system are considered as the achievement moments 
by the random process characterizing system, certain border, for example, zero. 
In [T5] are received approached exponential expressions for density of proba- 
bility of lifetime, accuracy of these expressions is estimated. In works [19l [20] 
lifetimes of molecules are investigated, the affinity of real distribution for life- 
time and approached exponential model is shown. It is possible to specify and 
other works (for example [5TJ [55J [53J ) where physical appendices of concept of 
lifetime widely applied in such mathematical disciplines, as reliability theory, 
the theory of queues and so forth (under names non-failure operation time, 
the employment period, etc.) are considered. In the present section lifetime 
joins in a circle of the general physical values, acting in an estimation or man- 
agement role (in terminology of the theory of the information [24]) for the 
quasi-equilibrium statistical operator that allows to receive the additional in- 
formation on system. In [53] it is noticed, that three disciplines grow together: 
statistical thermodynamics, Shennon's theory of the information and the the- 
ory of optimum statistical decisions. Accordingly, all correlations written down 
in the present work can be interpreted in terms of the theory of the information 
or the theory of optimum statistical decisions [55] . 

Let's notice, that in a case when value d\ng q (t — y, —y)/dy (the operator 
of entropy production a [1]) in the second item of the right part ([T]) does not 
depend from y and is taken out from under integral on y, this second item 
becomes cr(T), and expression (p} does not depend on form of function p q (y). 
There is it, for example, at g q (t) ~ exp{— at}, a = const. In work |26| ) such 
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distribution is received from a principle of a maximum of entropy at the set of 
average values of fluxes. 

3.2 Physical sense of distributions for past lifetime of sys- 
tem 

As is known (for example, [16], [18], [27]), exponential distribution for lifetime 

Pg{y) = eexp{-ey}, (3) 

used in Zubarev's works [TJ [5], is limiting distribution for lifetime, fair for 
large times. It is marked in works [1] [2] where necessity of use of large times 
connected with damping of nonphysical initial correlations. Thus, in works 
[TJ [2] the thermodynamic result limiting and universal is received, fair for all 
systems. It is true in a thermodynamic limit, for infinitely large systems. 
However real systems have the finite sizes. Therefore essential there is use of 
other, more exact distributions for lifetime. In this case the unambiguity of the 
description peculiar to a thermodynamic limit [28j is lost. 

For NSO with Zubarev's function ([3]) the value enter in second item 

-f p q {y)dy = exp{-ey} = l-ey+(eyf/2-... = 1-y / (T)+y 2 /2{T) 2 - ... (4) 

Obviously that to tend to infinity of average lifetime, (r) oo, correlation (01 
tends to unity. 

Besides exponential density of probability ([3]), as density of lifetime dis- 
tribution Erlang distributions (special or the general), gamma distributions 
etc. (see [HI [17]), and also the modifications considering subsequent com- 
posed asymptotic of the decomposition [57] can be used. General Erlang dis- 
tributions for n classes of ergotic states are fair for cases of phase transitions 
or bifurcations. For n — 2 general Erlang distribution looks like p q (y) — 
9pi exp{— piy} + (1 — 6)p2 exp{— P2y}, < 1< Gamma distributions describe 
the systems which evolution has some stages (number of these stages coincides 
with gamma distribution order). Considering real- life stages in non-equilibrium 
systems (chaotic, kinetic, hydrodynamic, diffusive and so forth), it is easy to 
agree, first, with necessity of use of gamma distributions of a kind 

p 9 (y) = e( £2 /)^ 1 exp{-e 2/ }/r(fc) (5) 

(X(k) is gamma function, at k = 1 we receive distribution l[3"])). and, secondly, 
- with their importance in the description of non-equilibrium properties. 

More detailed description p q (u) in comparison with limiting exponential 
allows to describe more in detail real stages of evolution of system (and 
also systems with small lifetimes) . Each from lifetime distributions has certain 
physical sense. In the theory of queues, for example |34] , to various disciplines of 
service there correspond various expressions for density of lifetime distribution. 
In the stochastic theory of storage [34], to these expressions there correspond 
various models of an exit and an input in system. 
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The value e without taking into account of dissipative effects can be defined, 
for example, from results of work [18]. The value e is defined also in work [3] 
through average values of operators of entropy and entropy production, flows 
of entropy and their combination. 

How was already marked, it is possible to specify very much, no less than 
1000, expressions for the distributions of past lifetime of the system. Certain 
physical sense is given to each of these distributions. To some class functions of 
distributions, apparently, some class of the physical systems corresponds, the 
laws of relaxation in which answer this class of functions of distributions for 
lifetime. 

3.3 Influence of the past on non-equilibrium properties 

A). Expressions for average fluxes. 

In [12] by consideration of the paradox connected with a waiting time, the 
following result is received: let X\ = Si; X2 = S% — Si;... are mutually in- 
dependent also it is equally exponential the distributed values with average 
1/e. Let t > is settled, but it is any. Element Xk, satisfying to condi- 
tion Sk-i < t < Sk, has density v t (x) — e 2 x exp{— ex}, < x < t;Vt(x) 
- [ I + ex) exp{— ex}, x > t. In Zubarev's NSO [TJ [2] the value of lifetime to a 
present moment t, belonging lifetime Xk, influence of the past on the present 
is considered. Therefore the value p q (u) should be chosen not in the form of 
exponential distribution ([3]), and in a form 

P q (y) = £ 2 yexp{sy}, (6) 

that in the form of gamma distribution |5j) at k = 2. In this case distribution 
(|6| coincides with special Erlang distribution of order 2 [16] , when refusal (in 
this case - the moment t) comes in the end of the second stage [16], the sys- 
tem past consists of two independent stages. Function of distribution is equal 
P q {x) = 1 — exp{— ex} — exexpj— ex},p q (u) = dP q {u) j 'du, unlike exponential 
distribution, when P q {x) = 1 — exp{— ex}. The behaviour of these two densities 
of distribution of a form © and © essentially differs in a zero vicinity. In case 
of © at system low probability to be lost at small values y, unlike exponential 
distribution ^ where this probability is maximal. Any system if has arisen, 
exists any minimal time, and it is reflected in distribution 

In work |10j it is shown, in what differences from Zubarev's distribution {1} 
with exponential distribution of lifetime ([3]) results gamma distribution ([5]), 
|6]) use. Additional items in NSO, in integral of collisions of the generalized 
kinetic equation, in expressions for average fluxes and self-diffusions coefficient 
are considered. The same in [10] is done and for special Erlang distribution k = 
2,3,4...,n,P 9 (x) = l-exp{-ex}[l + ex/l! + ... + (ex) fc -y(fc-l)!];e = k/(T). 
For distributions ([6]) is correct correlation (4), value — J p q (u)du — > 1 by 

Thus the multi-stage model of the past of system is introduced. Non- 
equilibrium processes usually proceed in some stages, each of which is char- 
acterized by the time scale. In distribution ^ the account of two stages, 
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possibly, their minimal possible number is made. Other distributions can de- 
scribe any other features of the past. Corresponding additives will be included 
into expressions for fluxes, integral of collisions, kinetic coefficients. Besides 
special Erlang distributions with whole and specified values k = n, that does 
not deduce us from set of one-parametrical distributions, the general already 
two-parametrical gamma distribution where the parameter k can accept any 
values can be used. In this case (r) = k/e. The situation (formally), when 
k < 1 is possible. Then sources will tends to infinity, as (t — to)| t Z^ to — * 00 
at k < 1. This divergence can be eliminated, having limited from below the 
value t — to of minimal lifetime value r ro ,- n , having replaced the under zero limit 
of integration on T m i n . Then to expression for a source |2]) it is added item 
[{eY min ) k - x /T{k))e exp{- £ r mm } In g q (t - 
B). Entropy production. 

Expressions for average entropy production received, for example, in [3], 
also depend from w (or p q {u) - in designations [3] and this work), i.e. on 
the chosen form of density of probability of distribution of the past of system. 
So, for average entropy production er = dS/dt in work [4] expression a(t) = 
J dt'Wit, t')(a(z\t, 0);a(z\t', t' — t)\t) k is received, where z are points in 
phase space , <j(z\t', t'—t) = —din g q (z\t' , t'—t)/dt', (<r(z\t, 0); a(z\t', t'—t)\t) k = 
(fc!)" 1 / dza{z\t, 0)a(z\t', t'-t) j£ dt x W{t, h)<T{z\t x ,H-t)... j£ dt k ^W{t 7 t k -x) 
cr(z\tk-i, ffe-i — t)Q q {z\t, 0),w(t, t') = dW(t,t')/dt' is "auxiliary weight func- 
tion" (in terminology [4]); w(t,t') it is designated above as p q (u);w(t,t') — 
p q {u = t — t'). For the limiting exponential distribution ([3j used in Zubarev's 
NSO, W(t, t') = exp{e(i' - t)}. 

4 Systems with infinite lifetime 

Above, as well as in work [TD], additives to NSO in the Zubarev's form are 
received for systems of the finite size, with finite lifetime. We will show, as for 
systems with infinite lifetime, for example, for systems of infinite volume, after 
thermodynamic limiting transition, the same effects, which essence in influence 
of the past of system, its histories, on its present non-equilibrium state are fair. 

In work |10j it is shown, as changes in function p q (u) influences on non- 
equilibrium descriptions of the system. But for those distributions p q (u), which 
are considered in [TU] (gamma-distributions, (HI, ©) the changes show up only 
for the systems of finite size with finite lifetimes. Additions to unit in equation 
([4| becomes vanishingly small to tend to infinity of sizes of the system and its 
average lifetime, as in the model distribution © used in Zubarev's NSO (JTJ) . 
For the systems of finite size and the exponential distribution results to nonzero 
additions in expression (fT]). Thus, these additions to NSO and proper additions 
to kinetic equations, kinetic coefficients and other non-equilibrium descriptions 
of the system, are an effect finitcness of sizes and lifetime of the system, not 
choice of distribution of lifetime of the system. We will find out, whether there 
arc distributions of lifetime of system for which and for the infinitely large 
systems with infinitely large lifetime an additional contribution to NSO differs 
from Zubarev'a NSO. 
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4a). Let's consider in quality p q {u) distribution of a form 

ku k ~ 1 p k 

Pq{u) = [i + m*] 8 ' 

received in work |17j , where k = 1 /r, f = — log p, t and v are parameters of 
scale and shift of logistical distribution f(x) = t _1 exp[(a; — is) / t] / {1 + exp[(x — 
v)/t]} 2 . In the correlatioon ([T]) the value J p q (u)du = — 1/[1 + appears. 
Average value of lifetime is equal 

up q {u)du = p" 1 B(l/(fc + 1), 1 - 1/fc), (7) 

(r) 2 = / u 2 p 9 (u)dM = /9 - 2 S(2/(fc + l),l-2/fc), 
Jo 

where B(, ) is beta function [32]. 

The value (r) in ([7]) to tend to infinity at a)p = 0, b)k = 1. Ratio of the 
second moment toward the square of the first is equal 

(T 2 ) fl(2/(fc + l),l-2/fc) 

<r> 2 s 2 (i/(fc + i),i-i/fc)' {) 

Expression © becomes vanishingly small with k — ► 1. Correlation (TTJ) for 
this distribution takes the form 

Ing(i) = / Pq(y)^Qq(t-y,-y)dy = lnQ q (t,0) + 
Jo 

1 dlnp 9 (i-u, -u) 
o + du )du - 

At k — ► 1 and finite values p we have a difference from the zero of additions to 
unit in expansion 



\ng(t) = lng q (t,0) + / (l-(up) k +(up) 2k -(up) 3k +...)(d\n g q {t-u, -u)/du)du. 
Jo 

But value k it is possible to define from correlation (JSj) , and, if ([8} is finite 
value not equal to the zero then k ^ 1. There is p — > 0, when additions becomes 
vanishingly small, as in the case of f4} of Zubarev's NSO. 

4b). Pareto distribution [T7] 



P 9 H = imw = T^-TTCTm' ^ = fc / fl - ( 9 ) 



/ca 70 k(k/p ) k 

[m+^F = [w+(fc/Po)P 

This distribution is received in [17j . as complex exponential distribution. It 
is supposed, that intensity p of exponential distribution f(u) = /oexp{— pu} 
represents random variable P with distribution fp(p). Then 

/•OO 

fr(u) =p q (u) = / pexp{-pu}f P {p)dp. 
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If to be set for function fp (p) by gamma distribution with density 

f P (p) = a k p k ~ l exp{-ap}/r(fc) (10) 
that we receive distribution ©, Pareto distribution. From © it is obtained: 



a , u 



k u 2 k{k + l) u 3 fc(fc + l)(fc + 2) 



P « {u)du = -J^T^ = - [1 -- ' 2a 2 6a 3 

(r)=a/(fc-l), k>l. 
We will consider two cases: 

a) . The parameter k is fixed, a = (T)(k — 1). Then, as in case of exponential 
((3|) or gamma distributions ([5]) for p q (u) additives to NSO are proportional to 
l/(r), they becomes vanishingly small with (r) — * oo as in 

b) . The parameter a is fixed. Then k = 1 + a/(r); — J p q (u)du = a k /(u + 
a) k = 1 - (u/o)(l + a/(T» + (u 2 /2a 2 )(l + a/(T»(2 + o/(T)) - (w 3 /6a 3 )(l + 
a/(r))(2 + a/(r))(3 + a/(r)) + ... -> 1/(1 + u/a), (r) -» oo,fc -► 1. We will 
mark that the second moment of Pareto distribution does not exist at k < 2 
P) . And at (r) -► oo,fc -> 1. 

Pareto distribution (J9]) corresponds to Tsallis distribution [30] 

1 1 _ 1 _ (fc + 2) _ (fc + 1) 

z[i + /3(«-i)«]V(?-i) ; fc+1 -(g-i) ; q ~(k + iy a ~ p ■ 

In Tsallis method averaging is conducted on distribution 

(fc/a)( fc+2 )/( fc + 1 ) 



P q {u) 



[1 + u/a 



fc+2 



Then 



(r) = Io l up T^ U = 7; tSt-tt-^; ^ = <r 2 >-<r> 2 ; 



J™ P i(u)du k' (r 2 ) (fc-i)' 



fc 



(s- 1)' (r 2 )' 7 " (u + a) fc 

_^ uk u 2 k(k + l) u 3 k(k + l)(fc + 2) 

a 2a 2 6a 3 

u m 2 2s m 3 2s(3s - 1) 

~ 0T + 2!(.s + i)(r} 2 ~ 3!(s + i) 2 (r) 3 + ■•' ' 

as in (|3|). By finite values of ratio s = D/{T 2 ) limiting behaviour of NSO is 
same, as in (0} and [TU] . 

4c). Let's consider one more distribution for p q (u), connected with degree 
laws. In work [12] by consideration of the renewal theory it is received distri- 
bution for length of an interval t — Sjy t , where S m are the renewal moments; 
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Sjv t <t< Sjf t +i- If to interpret the renewal moments as a birth and destruc- 
tion of system the interval t — Sw t represents time of past life of system, the 
value t — to from {I). In [T2] it is shown, that P{t — Sn ± > SW t +i — t > 
y} — > J x+ (1 — F(s))ds, where F(s) is distribution of random variable T, 

from sum = So + 2~i + ••• + 2n,/z = f °°{l — F(s))ds = f Q °° sF(s)ds. At 
large s: 1 — -F(s) ~ s~ a L(s), < a < 2, L(s) is slowly varying function. Thus, 
distribution p q (u) at large values t and a; + y looks like it~ a . We will break all 
time interval on two parts and we will describe at large times function p q {u) 
degree dependence, and on small times we set p q (u) gamma function of the 
form ©, ^TDJ) with k = 2, i.e ©. Thus, 

_ e 2 uexp{-ew}, u < c; , , 

where c is some value of time. From a normalization condition of distribution 
(fTTj) we find, that at 1 < a < 2, 

— (1 — a)eexp{— ec}(c+ 1/e) 

But value (T) disperses. Therefore we will be limited not to an infinite limit 
of integration on time, and some limiting value of lifetime r m . In this case the 
normalization condition gives value 

— (1 — a)eexp{— ec}(c+ 1/e) 

= (r m tt+1 - c-^ 1 ) ' 

and for average value of lifetime we receive expression 



2 lr (l-a) £ (c+l/ £ )(r m "+ 2 -c-"+ 2 ) (cV + 2c£+2) 

( )_ e +eXP| £C|[ (2-a)(T^+ 1 -^+i) £ J ' 

(12) 

If to fix parameters T m ,c,e, dependence e((r)), defined from (H2), will be 
positive for a limiting case interesting us 

T m -> oo, (r) -> oo, r m = lim(T}/T m 

at enough small values r m . For example, in case of extension exp{— ec} in series 
and restrictions of cubic items, we receive, that 

^ ((D-g) = (i-a)(r m "+ 2 -c-"+ 2 ) 

~ c 2 (c/3- 5 /2)' 9 (2-a)(r m Q+1 -c-«+ 1 )' 
In a limiting case 

r /-r\ 2 2 (.91 - r m) (1 - oQ ^ n 

F m -»oo, (r)->oo, e = , 31 = 7 > 0. 

&g x {a - 2) 

That was e 2 > 0, should be r m less g\. The ratio r m is finite at T TO ~ u~ a+2 as 
the value (T) disperses as u~ a+2 at u — > oo. Thus in a limiting case the value e 
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remains finite, and all additives entering in p q (u) and in additional expressions 
to NSO all additives are finite too unlike ((4]). 

If to fix parameters T m , c, e, defining dependence a from (r) from l|12p it is 
received, that at finite values c and e, at T m — > oo, (r) — > oo, r m = Hm(r)/r m 

(m + 2r m ) 

a = — — , to = exp^ — ec\[£c + 1 . 

(m + r m ) 

It is possible to consider and other limiting cases, and other distributions for 
Pq(u). But the general conclusion consists in that, as for infinitely large systems 
and infinitely large lifetimes the task of realistic distributions for time the 
lived of system a life changes a non-equilibrium state of system. The account 
of character of history of system, features of its behaviour in the past, can 
make essential influence on non-equilibrium behaviour of system in present 
time situation. 

4d). One more distribution for p q {u) can be received from results of works 
[31], integrating distribution P(E,T) on E 



p q {u) = exp{— 7u}(l — cexp{— 71*}) 1 ^ q ^; c = (q — l)jaq 1 ; (13) 



a = (expjjlPV} - l)" 1 ; (3=l/kT, 

where V = A is volume of metastable area, P is pressure, T is temperature 
|31j . The normalization of distribution (|13[) and its moments are expressed 
through incomplete beta- function [35]. For example, average value of lifetime 
is equal 



/n = a- 1 T 2 {\la) z F2{a- x ,a- 1 , l/(g -!);! + a" 1 , 1 + a~\ c) 
7<M ' 2 F 1 {a-\l/{q-l)-l + a-\c) 

^v i )(i + # F +^S F +...) 

y 1 + g(l+a) + 2q(l+2a) + —J 

r(a _1 ) is gamma function, n F m is hypergeometrical function [32]. The ratio 
of a dispersion to a square of average value is equal 

D/{T) 2 = 

WMb b b (^i) j 1 + b 1 + b 1 + b V(g + Sj c ) 

If to assume a little values 7 and 7a, and to be limited linear items, 
a za a 13 1 (1/a) 
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7<* 2 (r) a 2 7 _ 1 (i + g)(^-i) 

r 2 (l/a) ~ {l + afq a 

From here 7 — > at (r) — > 00 at finite values a, and we receive Zubarev's 
result ID, ([3]), ((I]) and [10], when additives to NSO becomes vanishingly small 
with (r) — y 00. The same gives also square-law approach on 7a. 

Let's consider now lifetime distributions of a various form to various time 
scales as (jlll) . In work [10 it was marked that during evolution the system 
passes various stages (kinetic, hydrodynamic, etc.). Lifetime can end at any 
stage. At different stages the functions g q accept a various kind. Therefore 
and expression for NSO (p} becomes complicated. 

4e). Let's consider distribution of kind 

. eexp{— eu}, u < c; 
P q (u) = { b ka» qt>r (14) 

(u+a) fc + 1 ' ' 

combining exponential distribution ([3]) for small times and fractional Pareto dis- 
tribution (|9|) for u > c. From a normalization condition 1 = e J Q exp{— eu}du + 
b f c °° a ( 1+ t/a) fe + 1 ^ ^ S f° un( i a normalizing constant, b = c k exp{—ec}, and then 
the first and second moments of lifetime: 

(T) = -[l-e- £C {l + ec)}~e- £C a[l+ '''' 



e L (1-fc) 

< r 2 > = ^ - e ^ + ec + — » ~ e a [(2^ - (T-fc) ~ 1] - 

For Pareto distribution at k < 2 the second moment does not exist [29] . 
The value (r) — > 00 with fc — > 1 from above. Limiting transition fc — ► 1 should 
be spent after thermodynamic limiting transition. The values e, c, a remains 
finite. Then 



exp{— eu}, u < c; 

: cxp{— ec} 
(l+u/a) k ' 



p q {u)du={ c " cxp{ - £c} 



p 9W cm t cexp {~ec}(l-u/a + ...), u > c. 

Thus, in this case additives to unit and to Zubarev's NSO are not equal 
to zero and for infinitely large systems with (r) — » cxd. Uncertain there are 
values of parameters e, c, a. As function p q (u) it is possible to choose both 
more simple and more difficult functions. 

4f). If to choose 

, , r eexpj— eu}, u < c; 
Hy ' 1 b, u> c, 

that b = °*y^~-c) ' integral in limits from to r m , but not from to 00. 
(T> = i[l-e-«(l + e c)] + e-"S^t^; 
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(r) -> oo, r r 



' r 2 ' 
(r 2 ) = | [i - e- £C (i + ec + ^)] + ^(r^ + r mC + c 2 ) ; 

(r 2 ) _ 2<T> 

r2 o - op ' m ' 

L rn ° m 

If after thermodynamic limiting transition T m — > oo, then — > 0, 

r eexp{— £it}, u < c; 

P,(W) = { cxp{- £C } > 

(r m - c ) ' " - c ' 

, x r exp{— £u}, u < c; 
H v ' 1 0, u > c. 

The value £ is hnite, and additives to unit are not equal to zero. But in this 
case the effect of "finite memory", limited on time by the size c, is observed. 
4g). For 

PiW-X be 2 exp{-e 2 u}, u>c, 

we write down from a normalization condition b = exp{— £ic}(l + e\c) and the 
two first moment 



1 2 (ftc) 2 

(r> = -e- £lC (l + £ lC )e- £2C (l + £ 2 c) + -[1 - e- £lC (l + £ lC + ^p-)]; 
£2 £1 I 

(T> -» 00 by £ 2 -> 0; 

<r»>=°[l-e— (l + ElC+ M! + ^] + 
£j 2 

+e- £lC 4e~ £2C (l + £ic)(l + £ 2 c + ^£L). 



Then 



/ \ j _ r (1 +ei")e U < c; 

P?W a I (l+£ic)e- £lC e- £2tl , u>c; 

p„(u)d M ^ . { ) ~ + -' u < c; 

F9W £2^0, (r) 00 V (1 + £ lC )e^ eiC , u>c. 

Additives to unit to tend to infinity of average lifetime are not equal to zero. 
4h). For 

£ 2 uexp{— eu}, u < c; 

[l+( 9 -l)a7exp{-7an}/ 9 ] 1 /(<!~i) I U - C ' 



combination of (J6J) and (|13jl . 
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7 a g - 1 7 a/g 1/a 1 ~ exp{-ec}(l + ec 
S(i_p,i)(l - l/(g-l),l/a) 

<r> = - [1 - exp{-ec}(l + ec + { -^-)]+ 

+ [1 - exp { -e C }(l + ^^(^ 

(a 2 j)2Fi(l/a, l/{q - 1); 1 + l/a;p) 

Bn_p ,i)(l — 1/(<Z — 1)? l/ a ) is incomplete beta function; (r) — > oo when 7 — > 0, 
as in 4d), (|13[) . 6^0 when 7 — > 0, after of thermodynamic limiting transition, 
and 

,-.,=> r e 2 uexp{ — eu}, u < c; 
H v ' 7—^0 1 0, u > c. 

, x , =>■ r eexpj— eu}(it + 1/e), u < c: 
p q {u)du { ^ 

y v ' 7^0 1 0, u > c. 

Thus, as in 4f) it is received "finite memory" on an interval (0, c), but contri- 
butions in NSO, amendments to Zubarev's NSO are finite and at (r) — > 00. 
"Finite memory" is possible and for "usual" functions of distribution by u > c 
and for fractional distributions of type of Pareto distribution. So, for 
4i). 

eexp{-eu}, u < c; 

Pq (u)du = { j,cxp{- 7 n} 

(1+au)" 1 ' " - C ' 

from a normalization condition 

a exp{— ec} exp{— 7a} 



6 = 



( 7 / a )i-mr(l - m,7/a) : 
F(, ) is incomplete gamma function, 



l r , , , r "fT(2 — m, 7/0) 

-[1 - exp{-ec}(l + ec + exp{-e C }[4^ '-UyL 

e a L (1 — m, 7/a) 



The first moment (r) — > 00 at a — > 0. In this case b — ► 7, and 

, . , =^> r eexp{— eit}, it < c; 
P<jl u J w a _^ q i ex p{— ec} exp{— 7u}, u>c. 

_ exp{-ew}, u < c; 

1 -7exp{-e C }/2^^, U >c; 

- / V (u)du { ex p{- £M l' u < c ; 

' 9 a^0 1 exp{— 7U} exp{— ec}, 14 > c, 



as in a case (|14[) . i.e. additives do not address in a zero at u > c. The value 
(r) — > co and at 7 — > co. In this case p q (u) = at u > c, memory is finite. 
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5 The conclusion 



As it is specified in work [8], existence of time scales and a stream of the 
information from slow degrees of freedom to fast create irreversibility of the 
macroscopical description. The information continuously passes from slow to 
fast degrees of freedom. This stream of the information leads to irreversibil- 
ity. The information thus is not lost, and passes in the form inaccessible to 
research on Markovian level of the description. For example, for the rarefied 
gas the information is transferred from one-partial observables to multipartial 
correlations. In work [3] values e = l/(r) and p q (u) = eexp{— eu} are ex- 
pressed through the operator of entropy production and, according to results 
[5], - through a stream of the information from relevant to irrelevant degrees 
of freedom. Introduction in NSO function p q (it) corresponds to specification of 
the description by means of the effective account of communication with irrel- 
evant degrees of freedom. In the present work it is shown, how it is possible 
to spend specification the description of effects of memory within the limits 
of method NSO, more detailed account of influence on evolution of system of 
quickly varying variables through the specified and expanded kind of density 
of function of distribution of time the lived system of a life. 

In many physical problems finitcness of lifetime can be neglected. Then 
e ~ l/(r) — > 0. For example, for a case of evaporation of drops of a liquid 
it is possible to show [33], that non-equilibrium characteristics depend from 
exp{y 2 }\y — e/(2A2) 1 ^ 2 , A2 is the second moment of correlation function of 
the fluxes averaged on quasi-equilibrium distribution. Estimations show, what 
even at the minimum values of lifetime of drops (generally - finite size) and the 
maximum values size y = e/(2A 2 ) 1 / 2 < 10~ 5 . Therefore finiteness of values (r) 
and e does not influence on behaviour of system and it is possible to consider 
e = 0. However in some situations it is necessary to consider finiteness of 
lifetime (r) and values e > 0. For example, for nanodrops already it is necessary 
to consider effect of finiteness of their lifetime. For lifetime of neutrons in a 
nuclear reactor in work [3] the equation for e — 1/ (r) which decision leads to 
expression for average lifetime of neutrons which coincides with the so-called 
period of a reactor is received. In work [35] account of finiteness of lifetime of 
neutrons result to correct distribution of neutrons energy. 

Use of distributions ([5]), J6]), ([9]), (fl3|) and several more obvious forms of 
lifetime distribution in quality p q (u) leads to a conclusion, that the deviation 
received by means of these distributions values \ng(t) from lng zu b(t) is no 
more e ~ l/(r). Therefore in expression ([1]) additives to the result received by 
Zubarev, are proportional e. This result corresponds to mathematical results of 
the theory asymptotical phase integration of complex systems |27| according to 
which distribution of lifetime looks \\kep q (u) = exp{— eu} + \ipi(u) + X 2 ip2(u) + 
... , where the parameter of smallness A in our case corresponds to value e ~ 
1 / (r) . Generally the parameter A can be any. 

For distributions of kind (lll|l , (fT4|) , (fTS"]) , having a various form for different 
times, additives to Zubarev's NSO are distinct from zero and for infinitely 
large systems with infinitely large lifetimes. For some distributions the effect 
of "finite memory" when only the limited interval of the past influences on 
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behaviour of system is observed. 

Probably, similar results will appear useful, for example, in researches of 
small systems. All greater value is acquired by importance of description of 
the systems in mesoscopical scales. A number of the results following from 
interpretation of NSO and p q {u) as density of lifetime distribution of system [3J , 
it is possible to receive from the stochastic theory of storage pH] and theories 
of queues. For example, in [34] the general result that the random variable 
of the period of employment (lifetime) has absolutely continuous distribution 
p q (u) = g(u,x) — xk(u — x,u),u > x > is resulted; g[u,x) — in other 
cases, where k(x,t) is absolutely continuous distribution for value X(t) - input 
to system. 

The form of distribution chosen by Zubarev for lifetime represents limiting 
distribution. The choice of lifetime distribution in NSO is connected with the 
account of influence of the past of system, its physical features, for the present 
moment, for example, with the account only age of system, as in Zubarev's NSO 
[TJ [U [3J at e > 0, or with more detailed characteristic of the past evolution of 
system. The received results are essential in cases when it is impossible to 
neglect effects of memory when memory time there is not little. The analysis 
of time scales as it is noted in [5], it is necessary to spend in each problem. 

Generalization of the received results on wider (generally any distributions 
F(x)) classes of distributions is received in work [T^] with use of methods of 
the renewal theory. In [12j it is shown, that the normalizing random variable 
(t — to)/t at t — » oo has limiting density g a (x) — (sinna / 'Tt)x~ a (1 — x) Q_1 , < 
a < 1, x € [0, 1], connected with functions of distribution F(x), having correctly 
varying tails, 1 — F(x) = x~ a L(x), < a < 1, where L(tx)/L(t) — > 1 at t — » oo. 
Average value (T)/t = (a — l)simra/sinir(a — 1). As (T)/t = 8 is small size 
values a are close to unit and 5 ra sin7ra/7r. At a w 1 — 5, sinil — S)tt = 
sinirS = Sir — (Sir) 3 /'3 + ... ~ ttS, we receive identity. At a w 1 — S distribution 
9a{x) ~ 6(1 — x)~ s /x 1 ^ 6 behaves in the similar image with eexp{— ex} at e ~ <5, 
differing at x — > 0. In this case universal distribution also is characterized 
only by one parameter a, but the limiting situation t — > oo and influence of 
tails of distribution, probably, not absolutely full describes past influence on 
the present as the near moments of time are thus more significant, with the 
memory which has not gone yet. 

In Prigogine's work [6] of function of distribution, evolving in course of 
time in accordance with the laws of mechanics, through transformation a dis- 
tribution function is put in accordance, the evolution of which is described by 
probabilistic rule. The role of such transformation in the method of NSO plays 
averaging on the density of distribution of time by the spent system of life. 

If type of source in Liouville equation for a non-equilibrium statistical oper- 
ator in the form of Zubarev [5] it is possible to confront with a linear relaxation 
source in Boltzmann equation, more difficult types of sources, got from other 
distributions for lifetime of the system, it is possible to compare to more re- 
alistic type of integral of collisions, that is explained by the openness of the 
system, by its co-operation with surroundings and finitcncss of lifetime of the 
system, and also coarsening for physically infinitely small volumes. 

The problems of correlation of spatial and temporal descriptions are inter- 
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csting. So, spatial description in works of Klimontovich [T3], smoothing out on 
physically infinitely small volume, yields to results similar with the results of 
Zubarev, to the same expression for a source. Account of influence of surround- 
ings, interaction with other systems (spatial task) also brings to the similar 
results. So, in works of MacLennan, short exposition of which and accordance 
with Zubarev's works is given in the Appendix 2 to the book of Zubarev [2J, 
a source in Liouville equation become formed by thermodynamical variables 
- temperature, chemical potential and speed characterizing surroundings, not 
details of his microscopic state. If, as in works [H[S] to conduct replacement of 
temporal argument of thermodynamical variables and to increase them on the 
proper "weigth functions" (which in [T] is interpreted as densities of distribu- 
tion of time by the spent system of life), we will get stated in the offered work 
results. 

Correlations between the exponential damping (|3|) and nonexponential func- 
tions just for small, non-Markov time scale factors, considered for a quantum 
case in works |36j . 
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